Stretching Semiflexible Polymer Chains: Evidence for the Importance of Excluded 

Volume Effects from Monte Carlo Simulation 
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Semiflexible macromolecules in dilute solution under very good solvent conditions are modeled 
by self-avoiding walks on the simple cubic lattice (d = 3 dimensions) and square lattice (d — 2 
dimensions), varying chain stiffness by an energy penalty e^ for chain bending. In the absence of 
excluded volume interactions, the persistence length t v of the polymers would then simply be l v = 
£b(2d—2)~ 1 q^ 1 with qt = exp(— et/ksT), the bond length lb being the lattice spacing, and ksT is the 
thermal energy. Using Monte Carlo simulations applying the pruned-enriched Rosenbluth method 
(PERM), both qb and the chain length N are varied over a wide range (0.005 < < 1, N < 50000), 
and also a stretching force / is applied to one chain end (fixing the other end at the origin). In the 
absence of this force, in d = 2 a single crossover from rod-like behavior (for contour lengths less than 
i p ) to swollen coils occurs, invalidating the Kratky-Porod model, while in d = 3 a double crossover 
occurs, from rods to Gaussian coils (as implied by the Kratky-Porod model) and then to coils that 
are swollen due to the excluded volume interaction. If the stretching force is applied, excluded 
volume interactions matter for the force versus extension relation irrespective of chain stiffness in 
d = 2, while theories based on the Kratky-Porod model are found to work in d — 3 for stiff chains 
in an intermediate regime of chain extensions. While for qb <C 1 in this model a persistence length 
can be estimated from the initial decay of bond-orientational correlations, it is argued that this is 
not possible for more complex wormlike chains (e.g. bottle-brush polymers). Consequences for the 
proper interpretation of experiments are briefly discussed. 



PACS numbers: 



I. INTRODUCTION 



The response of macromolecules with linear chemical 
architecture to mechanical forces pulling at their ends has 
been a longstanding problem in the statistical mechan- 
ics of polymers [l|-[29|. Particular interest in this prob- 
lem is due to advances in experimental techniques of sin- 
gle molecule measurements, probing the tension- induced 
stretc hing of biological macromolecules such as DNA [3(| , 
RNA 31], proteins [32[ and polysaccharides [33|. But 
also insight into the structure-property relationships of 
synthetic polymers, e.g. bottle brushes [Hj], has been 
gained by such experiments. However, despite extensive 
work on these problems, important aspects are still not 
well understood, even for the relatively simple case of 
macromolecules in dilute solutions of good solvent qual- 
ity, disregarding the interesting interplay of chain stretch- 
ing and collapse that occurs in poor solvents [2(| [35l - l38| , 
and also the inter play of chain stretching and adsorption 
on substrates |39l - [43 |. 

An important aspect of these problems is local chain 
stiffness. Traditionally, chain stiffness is characterized by 
"the" persistence length P, [§] but evidence has been pre- 
sented (43l - l45| that the traditional definitions are not use- 
ful under good solvent conditions, where excluded volume 
interactions create long range correlations with respect to 
the conformational properties of a macromolecule [H, [4(| . 
For stretched flexible polymers under a force / the stan- 



dard theory uses the concept of "Pincus blobs" of 
size £p = ksT/f, ksT being the thermal energy, pre- 
dicting a crossover from a Hookean regime, where the 
extension (X) for a force (applied in x-direction) scales 
as (X) = (R 2 )(f / dksT) in d dimensions, to a nonlin- 
ear power law (X) cx f 1 ^^ 1 (y being the Flory ex- 
ponent, v rj 3/5 in d = 3 dimensions). In contrast, 
for stretched semiflexible polymers excluded volume is 
widely neglected in the literature [l2|, EH, [l8l-|20l and us- 
ing the Kratky-Porod model (K-P model) [47], |48| simple 
analytic relations between force / and relative extension 
(X)/L, L b eing the contour length of the polymer, are 
derived 0, Hj. We recall that L = N b £ b where N b is 
the number of effective bonds of length l b connecting 
the effective monomeric units of the macromolecule, and 
if excluded volume interactions were absent, we would 
have, for N b — > 00 in the absence of the force /, the end- 
to-end distance of the polymer chains as (the index "0" 
refers to / = 0) 



(R 2 )o = i h L = 2£ P L = 



(1) 



l p being the length of a Kuhn segment, n 
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N b £b/£k being the number of such Kuhn segments form- 
ing the equivalent freely jointed chain P, Q. However, 
neither nor l v can be defined straightforwardly in the 
presence of excluded volume forces [43Tj46l . l49j . 

With recent large scale computer simulations, we have 
studied the combined effects of local "intrinsic" chain 
stiffness and excluded volume interactions on the confor- 
mational properties of polymers in the absence of stretch- 
ing forces, both for d = 3 [H, [44j and d = 2 [45| dimen- 
sions. The present study extends this work, giving a 
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detailed study of force-extension relations for both d = 2 
and d = 3, complementing our results also by investigat- 
ing fluctuations (X 2 ) — (X) 2 , (R"x) of the chain linear 
dimensions in the direction of the force and perpendic- 
ular to it. Whenever possible, a comparison with the- 
oretical predictions will be given. Since our study is 
based on modelling polymers as self-avoiding walks on 
square and simple cubic lattices, the main focus of our 
work is on the regime of low and intermediate forces (for 
very high forces, a more realistic description of the lo- 
cal structure and energetics of a polymer chains, such 
as bond length, bond angle and torsional potentials, be- 
comes important [HI). 

The outline of our paper is as follows: in Sec. II, we 
summarize the state of the art, discussing in particular 
the theoretical results we want to compare to. Sec. Ill 
briefly describes our model and the simulation technique, 
while Sec. IV reviews the properties of chains in the ab- 
sence of stretching forces. Sec. V describes the effects of 
stretching forces on conformational properties in d = 2 
and Sec. VI in d = 3 dimensions. Finally Sec. VII gives a 
summary and an outlook on experimental work, as well as 
on the related but more complicated problem of the "elec- 
trostatic persistence length" in polyelectrolytes (see 
for a review and further references). 

II. THEORETICAL BACKGROUND 

A. Force-Extension Curves for Flexible Chains 

Suppose we fix one end of an isolated polymer chain 
at the origin and apply a force / = (/, 0, 0) acting along 
the x-axis to the other chain end. This means, we add to 
the Hamiltonian of the chain a potential 

U = -fX (2) 

where X is the x-component of the end-to-end vector R 

-> Nb _ 

of the chain. Noting that R = ^2 Si, where a; = r - l+ i — r) 

i=i 

is the bond vector connecting monomers, at sites fj and 
Fj-fi, with \cii \ — t b the bond length which we assume as 
a constant, X can be rewritten as 

X = £ b J2cos^, (3) 

i=l 

where is the angle between at and the x-axis. For a 
model of freely jointed chains (FJC), a straightforward 
calculation of the partition function of the chain yields 
the force versus extension curve in terms of the Langevin 
function BBSHH (in d = 3 dimensions) 

(X) = e b N b £(f£ b /k B T) (FJC) , (4) 
C(0 = coth(C) - 1/C, ( = fl b /k B T. (5) 



Note that £(( < 1) ~ (/3 and hence one finds for small 
forces that 

(X) « X -l\N b flk B T = (R 2 ) f/(3k B T) , 

f£ b /(3k B T) (FJC), (6) 

where (R 2 )o = £ b N b for a freely joined chain. Eq. © 
would apply for Gaussian chains in the continuum for 
arbitrary large extensions, while for the freely jointed 
chain Eqs. dU, J5|) describe a saturation behavior for large 
/, when (X) approaches the contour length L = N b £ b , as 

{X)/L^l-k B T/fi b (FJC). (7) 

Of course, excluded volume interactions significantly 
modify the behavior described by Eqs. ((I])-©. Already 
in the absence of the force, Eq. (TTJ) is replaced by 

(R 2 ) = Ce b N 2v , N b -> oo (SAW) , (8) 

where C is a (non-universal, i.e. model - or system- 
dependent) constant of order unity, and the exponent 
v « 0.588 in d = 3 dimensions [52|, while v = 3/4 
in d = 2 dimensions @, [t| [53[ {Polymer chains behave 
like self- avoiding walks (SAWs)}. Treating the potential, 
Eq. as a small perturbation in linear response one 
can generally show that 

(X) = (R 2 ) f/(dk B T) for small/ (9) 

and hence the relative extension becomes in this regime 

{X)/L = C{L/t b f v -\fl b )/{dk B T) (SAW). (10) 

Comparing to Eq. ([5]) we note that the relative extension 
is enhanced by a factor C(L/4) 2ly_1 in comparison with 
the result for the freely jointed chain. 

While for the freely jointed chain the linear behavior, 
Eq. ([6]), smoothly crosses over to the saturation behavior 
(X)/L — > 1, Eq. (O, for the swollen coil there occurs an 
intermediate regime with a nonlinear relation between 
extension and force. This regime was first discussed by 
Pincus [5j in terms of the scaling ansatz 

(X) = {R 2 )l /2 F{{R 2 )l /2 /S, p ) (SAW) , (11) 

where £p is the radius of a "tensile blob" (also called now 
"Pincus blob"), 

£ P = k B T/f , (12) 

and F is a scaling function. Of course, this description 
makes only sense if 

4 « £p « (R 2 )l /2 - (13) 

since the scaling law for a blob (£p « £ b g v with g 
monomers per blob) breaks down when g is no longer 
very large; then a gradual crossover to the behavior of 
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a strongly stretched freely jointed chain must occur (ex- 
cluded volume then becomes irrelevant). For £p of order 

(R 2 )^ 2 , F( V ) behaves as (r, = (R 2 )l /2 /t P ) 



(14) 



so that Eq. (fTTj) reduces to the linear response results, 
Eq. ©. In the regime where Eq. (flUl) holds, the confor- 
mation of the chain is a stretched string of N b /g blobs, 
i.e. in order to obtain (A) cx L one must require that 
F(t~i) cx if~l v ~ x and hence 

(A) cx £ P N b /g cx (kBT/fib) 1 - 1 /^ (SAW) , (15) 

and hence one finds that in this regime the relative ex- 
tension varies as 



(X}/L <x (f£ b /k B T) 



l/u-l 



(SAW) 



(16) 



In a biopolymer context, this old result due to Pincus [|[ 
was recently "rediscovered" by Lam 54]. This scaling 
behavior can be made somewhat more explicit using the 
scaling description of the distribution P/v t (A) in which 
(in the absence of a stretching force f) a displacement X 

m 



between two end monomers occurs 



(X) = k B TdlnZ(f)/df (SAW) , 



(17) 



where the partition function Z(f) is (Zq is a normaliza- 
tion constant out of interest here) 



Z(f) = Z Q J d d RP Nb (X)eMfX/k B T) 



Using the ansatz [53 



P Nb (X) cx %) cx ^expI-DyVCi-")] . 



(18) 



(19) 



where y = X/(R 2 )l /2 , (j). = (1 - -y + ud - d/2)/(l - v), 
7 is a standard critical exponent [6J and D is a constant, 
Eqs. (IT71) . (HHJ) can be worked out numerically. 

Wittkop et al. [To| derived Eq. (ITB1) without explicit 
recourse to a blob picture. Wittkop et al. [l(| tried also 
to provide Monte Carlo evidence for Eq. (fl~6|) . both in 
d = 2 and d = 3, but they had to restrict their study to 
very short chains (20 < N b < 100). Morrison et al. HH 
argued that chain lengths of at least N b = 10 3 are neces- 
sary to provide clear simulation evidence for the Pincus 
tensile blob regime (described by Eas -JTB")) . dT^l) ). In 
fact, using N b = 6000 Pierleoni et al. [lj] succeeded to 
obtain evidence in favor of the Pincus theory Q for the 
chain structure factor under stretch. However, we are 
not aware of systematic tests of Eq. (fl~6| for very large 
Nb, as shall be presented here. 

A very interesting issue are also the longitudinal and 
transverse fluctuations, in the extensions of stretched 
chains. For freely jointed chains Titantah et al. [l6| de- 
rived 



(A 2 ) - {Xf 



N b f b [l-2£(C)/C-C 2 (C)} (FJC) (20) 



(21) 



which for large £ = fl b /ksT reduces to 

(A 2 ) — (A) 2 s» N b £ 2 C 2 = (k B T/f£b) 2 N b £ 2 (FJC) 

The transverse fluctuations becomes 

(Rl) = (Y 2 ) + (Z 2 ) 

= 2N b £ 2 b £(0/(~2N b £ 2 {k B T/f£ b ) (FJC) (22) 

where the last expression again refers to / — > oo. Of 
course, Eq. (|22[) differs substantially from a continuum 
Gaussian model of a chain (there the transverse linear 
dimensions are not affected by the pulling force at all). 

In the case where excluded volume is taken into ac- 
count, one obtains using again an approach based on 
Eqs. (fl~7l) - (fT9j) the approximate expressions [l6[ 

((A 2 )-(A) 2 )/(i? 2 ) = ^-(^) 2 + C- 2 (SAW), (23) 



si 



Sl 



(Rl)/(R 2 )v = C 2 Kc 2 / Sl - 1] (SAW) , (24) 
where the functions Si(C) & n d Cj(£) are defined as 



oo 

Si(C) = J dysmh((y)y 1+cl > exp[-Dy 



1/(1-^1 



1,3, 



(25) 



and 



d(0 = j dy C os\v{Cy)y l+ t'eM-Dy 1/(1 - v) ] , i = 1,2. 
o 

(26) 

Alternative approximate expressions were derived by 
Morrison et al. 25] using the self-consistent variational 
method due to Edwards and Singh [Hj]. 



B. Semiflexible chains in the absence of stretching 
forces: Chain linear dimensions and bond vector 
orientational correlations 

Following Winkler [2^, [56[ we first consider a chain 
with fixed bond length l b but successive bonds are cor- 
related with respect to their relative orientations, 



r 1 - 



{cii ■ a i+ i) 



t b t, t 



(27) 



being the angle between the orientation of two suc- 
cessive bond vectors. For this model, in the absence of 
excluded volume effects, the end-to-end distance is well- 
known 0, n 



, 9l „,,1 +t 2t t Nb - 1 

^ = N ^— + N- b — 2 



(28) 



In the limit N b — > oo the correlation function of bond 
vectors decays exponentially as a function of their chem- 
ical distance s, 

(ck ■ 3 l+s ) = £ 2 (cos8{s)) = £ 2 (cos6) s = £ 2 exp{-£ b s/£ p ) , 

(29) 
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where we have introduced the notion of the persistence 
length £ p [t| [5l| which becomes in this case 

4/4 = -ln«cos6l)). (30) 

In the case of semiflexible chains one has (cos 9) « 1 — 
(8 2 )/2, since (9 2 ) then is small, and hence one finds (for 
N b ->• oo ) du 

£ p ^2l b /(9 2 ), (i? 2 ) « 4AV 2 /(fl 2 ) =2£ p £ b N b , (31) 

so in this limit the Kuhn length 4 {Eq. (fT])} becomes 
4 = 2£ p , as was anticipated. 

When one considers now the limit 4 — > 0, A& — > oo, 
keeping L = 4-^6 as well as £ p finite, one obtains from 
Eq. (2gH [11] 

(i? 2 ) = 2£ P L{1 eM-L/£ P )}} , (32) 

which is nothing but the result that one could have de- 
rived directly from the Kratky-Porod model [U |H| for 
wormlike chains, 

L 

k f ,d 2 r,<, , s 

n= 2j { d^ )ds > (33) 
o 

where the polymer chain is described by the contour r(s) 
in continuous space. The bending stiffness k is related to 
the persistence length £ p as 

t 

n = £ p k B T, d=3, or n=^k B T, d = 2. (34) 

We note in this context the connection to the lattice mod- 
els that will be studied in the present work, where we 
study self-avoiding walks on square (d = 2) and simple 
cubic (d = 3) lattices, using a "penalty energy" e b if the 
chain makes a bend (by a 90° angle). Relaxing the ex- 
cluded volume constraint by considering a "non-reversal 
random walk" [571 ] . where only immediate reversals of a 
simple random walk model would be forbidden, we im- 
mediately conclude that 

(oos0} = l/[l + (2d-2))gb], q b = cxv(-e b /k B T) , (35) 

and hence one would obtain for q b — > from Eq. (j3Tj) 
that 

£ p /£ b = l/(2q b ) d = 2, or 4/4 = 1/(4%) d = 3. 

(36) 

At this point, it is interesting to recall that the present 
lattice model can be described by the Hamiltonian H = 
e bSi(l - o,% ■ o-il^b) = £ J2i( 1 - cos#j), plus excluded 
volume interaction, i.e. it is a discretized version of the 
Kratky-Porod model plus excluded volume, with angles 
Oi restricted to 9i = and 9i — ±90°, respectively. For a 
corresponding continuum model for large e small angles 
would dominate, however: putting 1 — cos 9i ~ 9f /2 the 



corresponding Hamiltonian would be H = ( £ b/2)X)j^f 
so one would conclude that k ( and hence £ p ) are simply 
proportional to e b . Eqs. (f35|) . ((36|) rather imply £ p oc 
exp(e/ksT); this effect is due to the fact that only large 
nonzero angles ±90° are permitted. 

Eq. (|32|) describes the crossover from the behavior of 
a rigid rod for L < £ p , where (R 2 )o — L 2 , to Gaussian 
chains for L > £ p , where Eq. ^ holds, {R 2 )o = 2£ p L. 
However, neither the exponential decay of the correlation 
function of bond vectors {Eq. (|29[) } nor the Gaussian 
behavior implied by Eq. (|32p remain valid when excluded 
volume effects are considered. 

On a qualitative level, insight into the effects of ex- 
cluded volume on semiflexible chains can be gotten by 
Flory-type free energy minimization arguments 0, 0, [HI, 
|59( . Consider a model where rods of length £k and diam- 
eter D are jointed together, such that the contour length 
L = N b £ b = n£k- Apart from prefactors of order unity, 
the second virial coefficient then can be estimated as (in 
d = 3 dimensions) 

v 2 = £ 2 D . (37) 

The free energy of a chain now contains two terms, the 
elastic energy and the energy due to the excluded volume 
interactions embodied in Eq. ([37|) . The elastic energy 
is taken as that of a free Gaussian chain, i.e. F e i 
R 2 /(£kL). The repulsive interactions are treated in mean 
field approximation, i.e. proportional to the square of the 
density n/R 3 and the volume R 3 . Thus 

AF/k B T a R 2 /(£ k L) + v 2 R 3 [(L/£ k )/R 3 } 2 (38) 

Minimizing A_F with respect to R, we find for L — > oo 
the standard Flory result 

R w (V4) 1/5 £ 3/5 = (4^) 1/5 (iV 6 4) 3/5 ■ (39) 

However, for finite L (or finite N b , respectively), Eq. (|3"9"]) 
applies only when the chain length N b exceeds the 
crossover length N£ or when R exceeds the associate ra- 
dius R*, 

N b >N b *, N: =£ 3 /(£ b D 2 ), R*=£ 2 k /D. (40) 

If N b < N b the contribution of the second term in 
Eq. (|38|) would still be negligible for R 2 £ k L, where 
AF/k B T is of order unity, and hence for Af b < A^* the 
first term in Eq. ([38| dominates, and hence the chain 
behaves like a Gaussian coil. However, this Gaussian 
regime only exists if N b ^> Njj° d — £ k /£b, the number of 
monomers per Kuhn length. For N b < A^' od , the chain 
resembles a rigid rod. Thus we predict (in d = 3) two 
subsequent crossovers: 

i?«L, AT b < Af h rod = 4/4 (rod -like chain), (41) 
R w (4i) 1/2 , A^' od <N b < N b (Gaussian coil) , (42) 
R w (4-D) 1/5 ^ 3/5 , N b > N* (R > R*) (SAW) . (43) 
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Of course, the intermediate Gaussian regime of Eq. (|42 



only exists if iV fc * » A£ od , or alternatively 



4 > D . 



(44) 



E.g., in the case of bottle brush polymers with flexible 
backbone chains and flexible side chains under good sol- 
vent conditions evidence has been presented for the fact 
that the stiffness of these wormlike chains is only due to 
their thickness [H, HU , and hence 4 and £ p are of the 
same order as D (disregarding the difficulty to define ei- 
ther 4 or £ p in this case properly) and thus the intermedi- 
ate Gaussian regime does not occur. On the other hand, 
for a large number of real semiflexible macromolecules 
under good solvent conditions the double crossover de- 
scribed by Eqs. (|41l) - (|4"3")l has been clearly observed [60| . 

On the other hand, the situation is completely different 
in d = 2 dimensions, where Eq. (|37l) is replaced by 



(45) 



since a rod of length Ik blocks an area of order l\ by occu- 
pation from a (differently oriented) second rod. Eq 
in d = 2 becomes [45| 



AF/k B T « R 2 /{i k L) + v 2 R 2 [(L/£ k )/R 2 } 2 . (46) 
Minimizing again AF with respect to R yields now 



fl«("2/4) 1/4 i 3/4 «4 /4i3/4 . d 



(47) 



where now the size L* = 4A fc * where Eq. l[4"7]) starts to 
hold is 



L* = 4, i.e. m = N, 



rod 



(48) 



Thus we note that a direct crossover occurs from rods 
to swollen coils (exhibiting statistical properties of two- 
dimensional self-avoiding walks), and no regime with in- 
termediate Gaussian behavior occurs. 

An important consequence of the excluded volume in- 
teraction is that they cause a much slower asymptotic 
decay of orientational correlations than described by 
Eq. (I29p occurs. For fully flexible chains one has a power 
law [46| 

(ck ■ a i+s ) ocs^, p = 2(1 - v) , 1 « s < N b . (49) 

Eq. (|4"9"|) has been verified by extensive Monte Carlo 
simulations of self-avoiding walks on simple cubic {d — 
3) [H,lill and square (d = 2) [45| lattices; we shall recall 
these results and extend them in Sec. IV below. 

For semiflexible chains we expect a crossover from ex- 
ponential decay {Eq. (|29|) } to the power law {Eq. (|5"Tj)} 
to occur near s = N£, i.e. we make the speculative as- 
sumption that 

(di ■ a i+s ) w l\ exp(-s4/4) , K s < JV b * , (50) 



s£h 

(2, ■ a i+s ) « exp(-A b *4/4)£ 2 (-^)-' 3 , AT* « a « N b . 

(51) 

Note that the prefactor of the power law in Eq. (|5T|) was 
chosen such that for s = N£ (where s£ b 
crossover to Eq. ([5Hf is possible. In d - 



= L*) a smooth 
2, where A b * = 

and hence 4 = £fe = 24, we note that the factor 



exp(-A h *4/4) 
finds that 



0.14 and using /3 = l/2inei = 2 one 



{Si ■ d l+s )/£ 2 » 0.14(24/4) 1 / 2 s- 1 / 2 , A b * « a « A fc 

and hence there occurs an increase of the amplitude 
of the power law with £ p . However, in d = 3, where 
N */ N rod = (2£ p /D) 2 the factor exp(-A*4/4) = 
exp(-2A*/A h rod ) = exp[-8(4/L>) 2 ] for large £ p will 
dominate and hence lead to a strong decrease of the am- 
plitude of the power law in Eq. (|5ip. Of course, the 
crossover at A fc * is not at all sharp but rather spread 
out over several decades in s, and hence the observabil- 
ity of Eqs. (f50"|) . (|5"Tj) is rather restricted. Note, how- 
ever, that ifli ■ di+ s ) always exhibits a single crossover 
only (for A fc * — > 00), near s = A fc *: while the radius ex- 
hibits two crossovers in d = 3 (at N b = Njj od and at 
N b = N£), Eq. (|50| does not exhibit any change in be- 
havior when s ps A™ d = 2£ p /£ b . Thus we predict that 
for very stiff and thin chains (for which £ P /D and hence 
A fc */A™ d are large numbers) one can follow the exponen- 
tial decay exp(— s4/4) over severa l decades. In d = 2, 



where — N™ , this is predicted to be impossible; 



rather one can follow the exponential decay only from 
unity to about 1/e. We shall discuss in Sec. IV a Monte 
Carlo test of these predictions. 



C. Stretching of semiflexible chains 



There exists a rich literature [12J, [13|, UMM |24H26] l28( 
where a force term, Eq. ([2]), is added to the Kratky-Porod 
Hamiltonian, Eq. (|33[) . to obtain 



H 



d 2 r(s) , 
ds 2 ■ 



l ds - f 



dx(s) 
ds 



ds . 



(53) 



We shall not dwell here on the exact numerical methods 
by which force versus extension curves can be derived 
from Eq. (|S"3"|) . but simply quote approximate interpola- 
tion formulas [HH (which are known to deviate from 
the numerically exact solutions at most by a few percent), 



L 4(1 



{X)/Lf 4' 



! = 3 (K-P model), 
(54) 

d = 2 (K-P model) . 

At this point, we remind the reader that k — £ p kgT in 
d = 3 while k = £ p k B T/2 in d = 2 {Eq. (|53|)}. For small 



k R T 



3(A) 
4 L 



[1-{X)/Lf 8 
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/, Eqs. (l54l) . (|55|) are compatible with the relations that 
one can derive by treating / via linear response, (X) = 
f(X 2 ) /k B T and hence ((X 2 ) = (R 2 ) /d = 2£ p L/d) 

li ' d = 2 ' 3 ( K - Pmodel )> ( 56 ) 
while for large / we find 

(X)/L »1- l/^4f£ p /k B T, d = 3 (K-Pmodel), 

(57) 

or 

(X)/L sa 1 - l/J8f£ p /k B T, d = 2 (K-P model). 

(58) 

A further quantity of interest is the "deflection 
length" (6^, [12, i-e. the correlation length of fluctua- 
tions along a semiflexible poly mer. In the presence of a 
strong force it is given by [19j 

\ = (f/k B T£ p y 1 / 2 , or X/£ p = {f£ p /k B T)- 1 ' 2 . (59) 

When ft p exceeds k B T, A hence becomes smaller than 
£ p , and in this limit one expects that excluded vol- 
ume indeed becomes negligible. However, when A be- 
comes of the order of the bond length £ b , it is clear that 
the continuum description in terms of Eq. (|53|) breaks 
down, the discreteness of the chain molecule becomes 
relevant [l9|, [22], [25[, and a crossover to the behavior 
of a freely jointed polymer occurs, as was described 
by Eq. ([7J. Toan and Thirumalai [28[ have empha- 
sized that all polymers under sufficiently high stretching 
forces should show a crossover from a force law of the 
type of the Kratky-Porod model, f - (X)/L cx f- 1 / 2 
{Eqs. ([57]) . (|55|)}. to the law of the freely jointed model, 
I - (X)/L oc f- 1 {Eq. ©}, and they argued that the 
crossover force between both descriptions is obtained by 
putting A = lb, providing evidence for this concept both 
by an analysis of experimental data and by a study of 
various models. However, here we are mostly interested 
in the regime where (X) / L is significantly smaller than 
unity, and excluded volume effects are still relevant. 

In particular, for weak forces we can combine Eq. © 
with the proper relations for the linear dimensions of the 
semiflexible chains, as discussed in Eq. (|39|) for d = 3 
and Eq. ([37} for d = 2, respectively. Thus, Eq. (1551) gets 
replaced by 

(X)/L oc (f£ P /k B T)(D/£ p ) 2 ^(L/£ p )^ 5 , L > L* = l b N£ 

(60) 

in d = 3, while in d = 2 we have 

(X)/Loc(f£ p /k B T)(L/£ p ) 1 / 2 , L>£ p , (61) 

where factors of order unity have been disregarded 
throughout. We note that in this regime the relations 
(X)/L versus f£ p /k B T vary much more steeply than pre- 
dicted by Eq. (|56l) if L/£ p is very large. Furthermore we 
note from Eqs. ([9l- (fT5)l that the nonlinear Pincus force- 
extension relation, Eq. (fTB"]). sets in when the size of the 



Pincus blob is of the same order as the coil size \J (R 2 )o ( 
Eq. ((391) for d = 3, or Eq. (@7J) for d = 2, in the absence of 
a force), defining a crossover length tjp c and associated 
force f c , 

{ P , c = ^f = (£ p D)^L^, d = 3, 

Jc 

or 

Cp, c - ^ /4 L 3/4 , d = 2, (62) 

i.e. for a force for which the extension (X) is of the same 
order as the coil size J (R 2 ) . As one expected for scaling 
theories, the crossover between the various regimes occur 
when all characteristic lengths \J (R 2 ), (X), £p c are of 
the same order. In the non-linear regime, according to 
the scaling ansatz, Eq. (TTTJ), Eq. ([To"]) gets replaced by 
(taking v = | in d = 3) 

(X) cx (fc B T//4r 2/3 (VD/^) 1/3 L, d = 3, (63) 
or (recall v = 3/4 in d = 2) 

(X)oc(k B T/fi b )- 1 / 3 (£ p /£ b )^ 3 L, d = 2. (64) 

Thus we see that in the nonlinear regime the persistence 
of the chains leads to an enhancement of the chain ex- 

1 /3 _____ 

tension by a factor £ p . Of course, the relation Eq. (jM)) 
can only hold if a Pincus blob contains many Kuhn seg- 
ments, i.e. now £p = k B T/f 3> £ p is required. This 
condition is nothing but the condition (X)/L < 1, in 
the case d = 2, as expected. In d = 3, however, we 
expect that the Pincus regime, as described by Eq. (|63[) . 
already ends when the size of a Pincus blob, £p = k B T/ /, 
equals the crossover radius R*, Eq. (|4U)) (remember that 
only for radii exceeding R* the excluded volume effects 
dominate) . For the crossover force 

/* = k B T/R* cx k B TD/£ 2 (65) 

we find from Eq. (|63l) that 

(X)/L cx D/£ p cx f*£ p /k B T . (66) 

Comparing Eq. ([66]) with Eqs. (j54l) . (l56l) . we see that 
indeed for / sw /* a smooth crossover from the Pincus 
behavior, as described by Eq. (1551 . to the Kratky-Porod 
law for wormlike chains for which excluded volume is 
negligible, can occur. We also note from Eq. (l65l) that 
in cases such as occur for bottle brush polymers |43l. Eij 
where chain stiffness is due to chain thickness, £ p oc D, 
we would have /* = k B T/£ p as in the two-dimensional 
case, and then the Pincus regime (which applies for 
k B T I \J (R 2 )o </</*) becomes much broader and eas- 
ier observable. 

Finally, we consider again the fluctuations in the ex- 
tensions of stretched chains (which were considered for 
flexible chains in Eqs. (|20p - (|2"6"]) already). However, the 
results known to us for semiflexible chains are some- 
what scarce (although the end-to-end distribution func- 
tion of the Kratky-Porod wormlike chain has been dis- 
cussed @, H3-H3)- Marko and Siggia QJ] obtained for 



(a) 



<X>/L 



1-hT/fh 

An 



(b) 



<X>/L 



i- w/fk 



d = 2 



e©Vrv 

^ ' i l-Jk B T/fl P chifins 
% P = h 



'freely 
jointed 



l ,°far 

.•r 



10" 



linear S a - <A' : > ; 
response I 

'Pincus blobs 



K _ 'D 



Kratky-Porod 
regime 



^=/\<R 2 > 



fh 



1 linear 
resp- 
onse 

10" 2 <ft 2 >„ / 



/ b 2 



/ I l-Jk B T/fh 
Kratky-Porod regime 



10"* 10" J 10"' 10" 



1 10 10 4 




freely 

nted 
chlains 



10"* 10"° 10"' 10 



1 10 10 4 



FIG. 1: Schematic plot of the relative chain extension, ({X)/L), versus the scaled force, (/£ p /fcsT), in a log-log scale for d = 2 
dimensions (a) and d = 3 dimensions (b). Broken vertical straight lines indicate various (smooth, not sharp!) crossovers in 
the response to the stretching force /. The first crossover occurs at very small forces, when the tensile length £ p = ksT/f 
becomes equal to the chain size y/ (R 2 )o in the absence of forces. In the first regime (to the left of this crossover) the extension 
is proportional to the force, (X) oc (R 2 )of /ksT (linear response regime). To the right of this crossover, the extension versus 
force curve follows the Pincus power law, (X) oc = / 1//3 (d = 2) or ~ f 2 ^ 3 (d = 3), respectively. In d = 2 this power law 

regime extends up to the crossover where the tensile length equals the persistence length £ p , while in d — 3 the power law ends 
already at an earlier crossover, £ p = R*, R* being the crossover radius where excluded volume statistics comes into play for 
semiflexible chains: then there exists a regime described by the Kratky-Porod model, {X)/L oc flp/ksT. For £ p smaller than 
l p , the extension approaches saturation according to the Kratky-Porod relation, 1 — \fkaT / fl p , while for still larger forces 
(when the deflection length becomes comparable to the bond length, a further crossover to the behavior expected for freely 
jointed chains ({X)/L w 1 — ksT/fii,)) occurs. 



the side-to-side excursions of the chain over a contour 
length s the result 



([?±(a)-?±(0)] a ) = 



exp[-s(f/nk B T) 



1/21 



/ 



(f/Kk B T)^ 



(67) 

where k is the coupling constant of the Kratky-Porod 
model {Eq. ([33])}. of course. For st b = L = N b £ b the 
result (R 2 ± ) = N b £ b (2k B T/f) is identical to the large 
force limit for the flexible chains, Eq. For small s 

Eq. (EZD yields ([r ± (s) - r x (0)] 2 ) = S 2 l\(k B T '/ 'k/) 1/2 • 

Given the fact that in the Pincus regime the picture 
of the chain conformation essentially is a stretched string 
of Pincus blobs (inside a blob excluded volume statis- 
tics prevails), we know that there occur of the order of 
n = N b (£ b /£ p )/g such blobs per string where g is the 
number of Kuhn steps (containing £ p /£ b monomers each) 
per blob. Remember that the Pincus blob has the radius 
£ p , and built as a self-avoiding walk of g steps of length £ p 
so that c> = i P g\ i.e. g = {£, P /t P ) 1,v = (k B T/ jl v fl v 
(here we disregard the factor 2 between the effective 
Kuhn step length £k and the persistence length £ p . Thus 
n = N b {£ b l£ p ){!£ P lk B Tfl v ) If this string of blobs 
would be completely stretched in a rod-like configura- 
tion, its lateral width would be simply ~ £p = 
£ p (k B T/ f£ p ) 2 . However, this estimate neglects the ran- 
dom statistical fluctuations that the string of blobs may 
exhibit in the transverse directions. We may consider 
the problem as a directed random walk where each step 



has a component £p in the +x-direction and a transverse 
component ±c£;p, where c is a constant (c <C 1). If we 
have n such steps (n — N b (g£ p /£ b )~ 1 ), we hence predict 
(Rl) = c 2 i\n = (?i v {k B Tlfl v f- x l v N b l h . Of course, 
this result can only apply if n is large enough so that 
cn 2 > 1, because cannot be smaller then £p, of 

course. Hence we would predict from these speculative 
scaling arguments 

(Rl) <xe P n = £ P £ b N b (k B T/f£ p ) 2 - 1 /\ n -> oo (68) 



<^i> = £ = e 2 p (k B T/f£ p ) 2 , small 



(69) 



In d = 2, this result should hold up to a force / = 
k B T/£ p , where = l v . Then Eq. (51 predicts (R 2 ^) cx 
£ p £ b N b = £ p L, i.e. there a smooth crossover to the re- 
sult oc k B TL/f derived from Eq. ((67)) occurs. In 
d = 3, however, Eq. (f6"8")) is supposed to hold only for 
£ P > R* cx £ 2 /D. 

At the end of this section, we summarize our findings 
for the force extension curves (Fig.[T]). The key to identify 
the various regimes is A COMPARISON OF LENGTHS, 
namely the "tensile length" £p = k B T/f needs to be 
compared to the various characteristic lengths of the un- 
perturbed chain. 

The simplest case actually occurs in d — 2 (Fig. [T]i). 

For £p > \J (R 2 )o we are in the regime of linear re- 
sponse, the extension (X) scales linearly with the force 
{Eq. (|6T])}. For £ P « y/(R 2 ) the extension (X) and 
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\J (R 2 )o are of the same order, linear response breaks 
down, and a (smooth!) crossover to the nonlinear Pin- 
cus regime occurs, (X)/L cx (/Ip/feT) 1 / 3 {Eq. (|64| |. 
The chain can be viewed as a stretched string of "Pincus 
blobs" of diameter £p (inside the blobs excluded volume 
statistics prevails). Near £p = £ p the extension (X) al- 
ready is no longer much smaller than the contour length 
L itself. Only the regime where the extension approaches 
its saturation value, from £p = £ p down to £p = 1%/lp, 
when the deflection length becomes equal to the bond 
length, the Kratky-Porod model holds {Eq. ([58])}. while 
for still larger forces (where the deflection length would 
be less than a bond length) the discreteness of the poly- 
mer chain causes a different relation between force and 
extension {Eq. ([7])}, as indicated in the figure. Of course, 
for flexible chains £ p — £b (actually £ p is completely ill- 
defined then) and the Kratky-Porod regime disappears 
altogether. 

For d = 3 dimensions the situation is more compli- 
cated, since for semiflexible chains without force an- 
other regime appears, for distances in between £ p and 
R* = lp/D, where Gaussian statistics prevails, and 
this regime finds its correspondence in the force ver- 
sus extension curve. Thus, for very long semiflexible 
thin chains there are three regimes, where the force- 
extension curve exhibits power laws: for very weak forces 
(£p > \J (R 2 )a) the linear response regime occurs with 
(X) cx (R 2 )af /k B T {Eq. (pD|) }. then a nonlinear regime 
with (X)/L cx (/£p/fc s T) 2 / 3 (€ p /i?*) 1 /3 {Eq. (jn3D> fol- 
lows for y/ (R 2 )o > C p > R* cx £ p /D, and then the linear 
regime as described by the Kratky-Porod model follows, 
(X)/L cx f£ p /k B T {Eq. (OH)}, for R* > £p > £ p . For 
£p < £p the approach of (X}/L to its saturation value 
unity proceeds in a similar manner as in d = 2 {Eqs. (1571) . 
©}• 

Of course, the description in Fig. Q] contains the force- 
extension curves of fully flexible polymers as a limiting 
case: there both R* and £ v tend to £b, and the regime 
where the Kratky-Porod model is applicable gradually 
disappears. The same statement applies to semiflexible 
chains in d = 3 when stiffness is due to chain thickness, 
so that R* tends towards £ p , and the nonlinear Pincus 
blob regime {{X)/L cx (f£p/k B T) 2 ' z takes over and holds 
down to £ p ss £ p . Conversely, if one considers not very 
long chains, such that (R 2 )o < R* 2 , the regime domi- 
nated by excluded volume effects disappears from the pic- 
ture, and the Kratky-Porod model description becomes 
valid down to arbitrarily small forces. 

We end this section with several caveats: (i) All 
crossovers in Fig. [T] are smooth and we do not expect any 
sharp kinks at the crossover values of £p that are indi- 
cated by the vertical broken lines; rather gradual changes 
will occur on the log-log plot (X)/L vs. f£ p /k B T, spread 
out over (at least) a decade in f£ p /k B T. Consequently, 
a/ (R 2 )o must exceed £ p by four (or more!) decades, 
in order to resolve the multiple crossovers of Fig. [T] in 
d = 3. (ii) We have disregarded any special structures 
of the polymer such as oc-helices known for proteins, 



double helix-portions of copolymers formed from double- 
stranded DNA and other biopolymers, etc. [68, 69]. Any 
such special structures of biopolymers will lead to non- 
universal special features of the force-extension curve, in 
particular in the regime of rather low forces, but these 
are outside of consideration here. Rather only a generic 
description of the universal behavior of very long flexible 
or semiflexible polymers is within our focus, (iii) Very 
long polymers exhibiting a (swollen) random coil config- 
uration are expected to contain knots (to precisely define 
them, one can transform the polymer configuration into 
a closed loop by adding the end-to-end vector as an ex- 
tra special bond) [70j . Pulling such a chain at both ends 
will have the effect that the knots tighten, and the knots 
can only be made to disappear by moving them to the 
chain ends. Effects due to knots |7il |72| clearly are be- 
yond the realm of our scaling description (while in the 
computer simulations presented in Sees. IV- VI knots are 
automatically included implicitly in our models, though 
we do not make any explicit attempt to study their ef- 
fects). In view of all these caveats, the extent to which 
the scaling theory sketched in Fig. 1 is practically useful 
is a nontrivial matter. 



III. MODEL AND SIMULATION TECHNIQUE 

The model that we study in this paper is the classical 
self-avoiding walk (SAW) [g, [tl H3] on square and simple 
cubic lattices, where bonds connect nearest neighbor sites 
on the lattice, and the excluded volume interaction is 
realized by the constraint that every lattice site can be 
taken only once by an effective monomer occupying that 
site. We take the lattice spacing as our unit of length, 
£b = 1. Variable chain stiffness (or flexibility) then is 
introduced into the model by an energy eb that occurs 
for any kink (that is at a right angle and costs eb). 

No energy arises for 9 = 0° , of course, and hence in the 
statistical weight of a SAW configuration on the lattice 
every kink will contribute a factor qb — exp(— e/k B T). 
In the presence of a force /, the potential U written in 
Eq. ([2]) yields another factor b x with b — exp(//fcpT) 
to the statistical weight (as in Sec. II, the force / is as- 
sumed to act in the positive x-direction, and X is the x- 
componcnt of the end-to-end vector R of the chain). So 
the partition function of a SAW with Nb bends (Nb + I 
effective monomers) and Ab cn d local bends by ±90° is 

ZN,N hend (<lb,b) = C(N b ,N bcnd ,X)q b N ^b x . (70) 

config 

We have carried out Monte Carlo simulations applying 
the pruned-enriched Rosenbluth Method ffMf^ (PERM 
algorithm) using chain lengths up to Nb = 50000 in d = 3 
and Nb = 25600 in d — 2, varying also the chain stiffness 
over a wide range (0.005 < qb < 1.0). As mentioned 
already in Sec. II {Eq. (l3"rJl) }. this means that the persis- 
tence length £ p varies over about two orders of magnitude 
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FIG. 2: Semi-log plot of (cos#(s)} versus the contour length s, for qb in the range from qb = 0.1 to g& = 1.0 (a) (c), and for 
rather stiff chains, 0.005 < q b < 0.05, (b) (d). Data are taken for N b = 25600 in d = 2 (a) (b) and for N b = 50000 in d = 3 
(c) (d). The straight lines indicate fits of the initial decay to Eq. (|29[) . (cos#(s)) oc e~ alb ^ lp \ for flexible chains (g;, = 1.0 and 
q b = 0.4) meaningful fits are not possible. Estimates for E p /£ b are listed in Table HI and ITu Note the difference in ordinate scales 
between d — 2 and d = 3: in d = 2, deviations from Eq. ((29} have set in when (cos#(s)) has decayed to 1/e, in d = 3, however, 
for very stiff chains one can follow the exponential decay for almost two decades. 



(note, however, that in the presence of excluded volume 
one has to be very careful with the notion of a persistence 
length, particularly in d = 2 dimensions p3l - |45| V 



IV. SEMIFLEXIBLE POLYMERS IN THE 
ABSENCE OF STRETCHING FORCES 

In this section, we summarize our Monte Carlo results 
for bond orientational correlations and chain linear di- 
mensions obtained for the model described in the previ- 
ous section. While some of these results have recently 
been described in our earlier work [43l - l45| . the informa- 
tion provided will be crucial for the understanding of our 
results for the extension versus force curves as well. 

We start with the bond orientational correlation func- 
tion (cos (9(s)), Figs.[2El since the decay of this function 
with s is traditionally used to extract "the" persistence 
length l p , using Eq. (l29|) . As expected from Sec. II B, 
however, one must not rely on Eq. (|2TJ)) to describe the 
asymptotic decay of (cos#(s)) (in the limit where first 
Nb — > oo has been taken, so that one can study large s 



without being affected by the finite size of the chain) for 
large s, but rather one must consider the initial decrease 
of (cos 8(s)} with s, cf. Eq. (JSOJ. Since s = 0, 1, 2, 3, . . . 
is a discrete variable, such a fit becomes ill-defined for 
flexible chains; then the only possible procedure is to 
use Eq. (I3T)]) as a definition of the persistence length, 
l Vt B = — lb I ln[(cos#(s = 1))]. Both estimates for l p 
(from an extended fit over a range of s, and from the 
latter formula) are collected in Table HI together with 
the prediction based on Eq. (|36[) . where excluded volume 
is neglected. One recognizes that Eq. (I31)|) becomes ac- 
curate for d = 3 as the chains become very stiff, qb — > 
while in d = 2 Eq. (|3"6) {predicting l v jlb = O.bq^ 1 in this 
case} never becomes valid. As a consequence, we empha- 
size that the rule (based on the Kratky-Porod model) 
that for the same bending stiffness n/ksT (the contin- 
uum analog of our parameter qb) the persistence length 
in d = 2 is twice as large as in d — 3 is not accurate 
(since this rule fails only by about 24 %, in experimen- 
tal work where it was tried to extract estimates of £ p 
from adsorbed semi-flexible chains on two-dimensional 
substrates this problem was not noticed, due to other 
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FIG. 3: Log-log plot of (cos 9(s)) versus s, for — 0.1, 0.2, 0.4 and 1.0, including only data for Nt ~ 25600 in d = 2 (a) and for 
Nt = 50000 in d — 3 (b). The straight line indicates a fit of the power law, Eq. (|49[) , to the data for qb = 1.0, including only 
data for sib > 10 in the fit, and requesting the theoretical exponent, f3 — 2 — 2l>, with ^ = 3/4 in d — 2 (a), and v = 0.588 in 
d = 3(b). 
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FIG. 4: Semi-log plot of (cos(9(s)) versus the scaled distance s* — slb/lp along the chain, for d = 2 (a) and d = 3 (b). Data for 
£p/£b extracted from Fig. [2] as described above and listed in Table HI and ITTl were used. 



uncertainties in the data analysis). 

Fig. [3] plots our data for the bond orientational corre- 
lations in a log-log form, to clearly demonstrate that the 
asymptotic decay is a power law {Eq. (I49l| rather than 
exponential {Eq. <J29])}. Eqs. ([5D|). ((52)) in fact suggest to 
study (Sj • Oi+ s ) not simply as a function of s but as a 
function of the rescaled variable s* = s£b/£p- In d = 2, 
one expects a data collapse on a universal master curve 
and this is indeed found (Fig. 2£a)). No such simple 
scaling is possible in d = 3, however, as expected from 
Eq. 

Another measure of a persistence length in our model 
is the average number (n str ) of successive bonds along the 
chain that have the same orientation without any kink. 
The distribution P(n s t T ) of such straight sequences along 
the chain is plotted in Fig. [SJ We find that irrespective 
of dimensionality and for all values of qb the distribution 
shows a simple exponential decay 

f(n st r) = a p exp(-7j str /n p ), (71) 



and either the average (n s tr) = X) -P^stiO^-str or the 

»str=l 

decay constant n p can be taken as a characteristic (sim- 
ilar but not identical to the persistence length) of local 
intrinsic chain stiffness. Clearly, fitting the data shown 
in Fig. [S] is less ambiguous than fitting the data for the 
bond orientational correlations. We also note that (n stI ) 
has an obvious physical correspondence in real macro- 
molecules: in alkane-type chains, where the torsional po- 
tential has one deeper minimum (the "trans" state, tor- 
sional angle tp = 0°) and two less deep minima (gauche 
±, tp = ±120°) separated from the trans state by high 
energy barriers, n str simply is the number of successive 
carbon-carbon bonds in an all-trans configuration. Of 
course, in this case successive bonds in this state are not 
oriented along the same direction, since the ground state 
configuration of the alkanes is a zigzag-configuration, and 
so Eqs. (|2T))) , (|50"|) need to be generalized (bond angles 
need to be measured relative to their values in the "all- 
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TABLE I: Various possible estimates for persistence lengths, l p /£b from Eq. (|29[) . £p,e/£b from Eq. I|30[) . n p and {n str } from 
Eq. (|71[) including the fitting parameter a p , and £ p ,n/£b from Fig. [6l for semiflexible chains in d — 2 with various values of q;,. 



(//. 



0.005 0.01 0.02 0.03 0.05 0.10 0.20 0.40 1.0 



t p ,e/h 

a p 

n p 
(rw) 



123.78 62.20 31.36 

118.22 59.44 30.02 

0.009 0.017 0.034 

116.34 58.81 29.78 

118.06 59.78 30.48 

3.34 2.37 1.70 



21.55 13.41 7.53 4.16 2.67 

20.21 12.35 6.46 3.50 2.00 1.06 

0.051 0.085 0.168 0.331 0.646 1.539 

20.07 12.28 6.44 3.50 2.01 1.08 

20.69 12.85 6.97 4.02 2.54 1.64 

1.39 1.09 0.81 0.61 0.48 0.39 



TABLE II: Various possible estimates for persistence lengths, £ p /£b from Eq. (|29|) . £ p ,e/lb from Eq. (|30[1 . n p and (n str ) from 
Eq. (|71|) including the fitting parameter a p , and £ p ,n/£b from Ref. [44J] . for semiflexible chains in d — 3 with various values of 



lb 



0.005 0.01 0.02 0.03 0.05 0.10 0.20 0.40 1.0 



£ P /h 
£ P ,e/£b 
a p 
n p 
(n st v) 

tp, R /£b 



52.61 26.87 13.93 9.54 5.96 3.35 2.05 - 

51.52 26.08 13.35 9.10 5.70 3.12 1.18 1.12 0.67 

0.02 0.04 0.08 0.12 0.19 0.38 0.73 1.42 3.37 

51.17 25.95 13.30 9.07 5.68 3.12 1.82 1.13 0.68 

51.72 26.50 13.83 9.60 6.20 3.65 2.36 1.70 1.29 

5.35 3.49 2.39 1.94 1.54 1.12 0.87 0.71 0.61 



trans" configuration). Similarly in biopolymers a se- 
quence of n str bonds in an ev-helix configuration can be 
the right object to characterize stiffness. Such consid- 
erations will be needed when one wants to adapt our 
findings to real polymer chains. Finding an analogue of 
Eq. (fTTj) that is generally valid for off-lattice models is in 
interesting problem but beyond the scope of the present 
study. 

We now turn to the end-to-end distance of the chains 
(Figs. |51 ©• As predicted by the theoretical considera- 
tions of Sec. II, we find in d = 2 and d = 3 dimensions 
very different behaviors. In d — 2 (Fig. |6]) the data for 
small enough Nb show the rod-like behavior, (R 2 ) oc , 
indicated by the slope of the straight line in the left of 
Fig-E^a), (b), and then (R 2 e ) oc with v = 3/4 reaches 
a broad maximum, and thereafter decreases only a little 
bit and then settles down at the limiting value expected 
for two-dimensional self-avoiding walks. As the rescaled 
plot (Fig. [Sib)) shows, there is a single crossover from 
rods to self-avoiding walks, and irrespective of stiffness 
there is never a regime where the Gaussian plateau pre- 
dicted by the Kratky-Porod model {Eq. ([32)) } describes 
part of the data approximately. Of course, the latter can 
describe the initial rod-like behavior 49] but this is of 
little interest and clearly from this regime one cannot es- 
timate £ p reliably at all. Interestingly we find from the 
rescaled plot (Fig.^b)) that the maximum which appear 
at Nb — A^ nax in Fig. M,&) rather accurately coincides 
with the value Nb = N£ , the chain length corresponding 
to the effective Kuhn length £k — 2£ p . As an immediate 
consequence of this finding we can suggest as a recipe 
for experimentalists who analyze end-to-end-distances of 
two-dimensional adsorbed chains to plot their data in 



analogy to Fig. [UJa)): if their chain lengths Nb are long 
enough to reach the region where the maximum jV" lax 
in such a plot occurs, they can immediately estimate the 
persistence length as 

£ P = £ b Nr x /2 ■ (72) 

For this method to work, it is not necessary at all to 
have chains long enough to see the asymptotic d = 2 SAW 
behavior, (R 2 ) oc N^ 2 . Since the Kratky-Porod model 
{Eq. (J22)} has so widely been used by experimentalists to 
fit their data and by theorists to build more sophisticated 
extensions on it, we emphasize again that Eq. Q32[) is ac- 
curate only in the rod-like regime and in the initial part 
of the crossover towards the self-avoiding walk regime as 
shown in Fig. [l^b). Note that unlike experimental work, 
this comparison does not involve any adjustable param- 
eter whatsoever. We see that for small Nb and small 
qb theory and simulation agree qualitatively, but in this 
regime, where curves for small Nb collapse on the straight 
line (R 2 }/£bNb — £bNb, and then gradually bend over to 
a slower increase, the data are not very sensitive to the 
actual value of £ p . For N b < 2iV b * the Kratky-Porod 
result slightly overestimates the actual data, while for 
Nb 2A^ it strongly underestimates them, since the in- 
crease proportional to (R 2 ) oc cannot be described. 
Clearly, the plateaus predicted by Eq. (J32]) for Nb > 2N£ , 
as displayed in Fig.^b), do not have any correspondence 
to the actual data. 

However, in the three-dimensional case the situation 
is clearly different [HI, IZij . We shall not reproduce in 
full detail the data published already elsewhere [43|, |44| 
but only show as a summary of the scaling plots in the 
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FIG. 5: Semi-log plot of the distribution P(n str ) of n str successive bonds along the chain which continue straight along a lattice 
direction until a kink appears, versus n str for rather flexible chains, i.e. 0.1 < qb < 1.0 in d = 2 (a) and d = 3 (c), as well as 
for rather stiff chains, i.e. 0.005 < qb < 0.05, in d = 2 (b) and d = 3 (d). The straight lines indicate fits to simple exponential 
functions, Eq. (|7ip. P(n str ) = a p exp(— n stl /n p ), with constants a p and n p quoted in Table U and [TTJ All data are taken for 
N b = 25600 in d = 2 and % = 50000 in d = 3. 



Kratky-Porod representation for the three-dimensional 
cases (Fig. [7]). Now there is clear evidence from the 
data (Fig. E{a)) that with increasing stiffness (decreas- 
ing qb) a Gaussian plateau in the plots of (R^) / (2lbNbl p ) 
versus N b /Nl od {q b ) develops, before the regime ruled by 
excluded volume interactions sets in. Here we rescale the 
chain length Nb such that data collapse occurs in the rod- 
like regime, i.e. Nb is rescaled with N b od = 2£ p /£ b - We 
see that with increasing l v the data gradually approach 
the Kratky-Porod result {Eq. (J22)} over an increasing 
range of Nb/Njj od , while ultimately the data increase be- 
yond the Kratky-Porod plateau values, to cross over to 
the asymptotic relation (Rl)/N b oc N^^ 1 as it should 
be {compare Eqs. Alternatively, we can es- 

timate another crossover chain length N b (q b ) such that 
the curves collapse in the regime of large Nb , so that the 
asymptotic regime where excluded volume interactions 
dominate, shows proper scaling behavior (Fig. [7p). Ob- 
viously, while in d = 2 dimensions N b (q b ) — N b od (q b ), 
so there is no need to distinguish these crossover chain 
lengths at all, (Fig. and there is a single crossover 
from rods to self-avoiding walks described by one univer- 
sal crossover scaling function, this is not true in d = 3 



dimensions: there occur two successive crossovers, from 
rods to Gaussian coils at Nb = N b ° d , and from Gaus- 
sian coils to three-dimensional self-avoiding walks, at 
Nb = Nj* . Of course, these crossovers are rather gradual 
and not sharp: therefore a well-defined Gaussian plateau 
comes into existence only for N b * 3> A r ™ d , which requires 
extremely stiff chains. These findings are in beautiful 
qualitative agreement with the theoretical considerations 
of Sec. II B and with available experiments in d = 3 
that did show two successive crossovers [60[ • Surprisingly, 
some authors [76| claim to have observed two successive 
crossovers (with an intermediate Gaussian regime) for 
two-dimensional adsorbed chains. We suspect that the 
observations may be due to incomplete equilibration of 
the chains, and we feel that the theoretical interpretation 
given there is inappropriate, however. 

It remains to test to what extent the predictions given 
in Sec. II B for the crossover chain lengths N™ d ,N b * are 
actually compatible with our data. First of all, Fig. Uta) 
illustrates that both in d = 2 and d = 3 the region where 
lb follows the simple asymptotic power law £ b oc q^ 1 is 
quickly reached, and we have ample data where £ p ex- 
ceeds lb by at least an order of magnitude. 
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FIG. 6: Log-log plot of the rescaled mean square end-to-end distance in d = 2 dimensions, (Rt) /(2£lN^ 1/ ) with v — 3/4, versus 
Nb (a) and the normalized rescaled mean square end-to-end distance (Re) / (2£ p £bNb) versus rescaled chain length Nb/Nl° d (b). 
These data include chain lengths Nb up to Nb = 25600, and all values of the stiffness parameter qb, as indicated. Straight lines 
in (a) show the slope 2 — 2v = 0.5 describing the rod-like regime (that occurs for small Nt) and the slope 1 — 2v = —0.5 that 
would occur if a Gaussian-like regime was present (which is not). Dotted horizontal plateaus for large Nb in (a) show estimates 
for £p,n(qb)/£b (Table H]). Part (b) shows that all data collapse to a single master curve which describes a crossover from a 
rod-like regime to a self-avoiding walk regime. The Kratky-Porod function, Eq. ((32}, indicated by the dotted curve (WLC) is 
also shown for comparison. The chain length N£ = N b od = £k/£ b = 2£ p /£ b describing the number of bonds per effective Kulm 
segment is extracted from the persistence length estimates (Table H}. 




FIG. 7: Log-log plot of {R 2 e )/(2£ b £ p (q b )N b ) versus N b /Nl° d {q b ) with Nl° d = 2£ p (q b )/£ b . As always, £ p (q b )/£ b is extracted from 
the initial decay of the bond vector correlation function (Fig. [2] and Table ITTj) . Now the initial part of the data scale, and the 
smaller q b (increasing £ p ) the more the data follow the Kratky-Porod function, Eq. (|32l) . indicated by the dotted curve labeled 
WLC. Part(b) shows the same data but plotted versus N b /N b (q b ) (a) where N b (q b ) is defined such that for large Nb an optimal 
data collapse on the straight line representing the three-dimensional self-avoiding walk behavior (Re) oc N 2v is obtained. Note 
that data that fall below the Kratky-Porod plateau (horizontal straight line) for different stiffness parameters systematically 
splay out, there is no scaling over the full parameter range. 



Fig. [HJd also illustrates that in the case of a sin- 
gle crossover Eq. (|T7|) is quantitatively verified, since 
Eq. g?D says R 2 oc 4 /2 L 3 / 2 = f b /2 el /2 N^ /2 , and hence 
using (R 2 ) = 2£ Pt u£bNb we would conclude £ p _r — 
(ipib) 1 / 2 /\/2, if the proportionally constant in Eq. (|47|) 
is taken to be unity. Of course, only the exponent in 
this relation and not the prefactor can be taken seriously. 
However, in d = 3 the theoretical relations N£ oc & and 

C-p,r oc £^ 5 {Eqs. (PffJl) . (02])} are not quantitatively ver- 



ified: rather we found effective exponents £ p ^r oc l®' 56 
and oc £ p 5 (Fig. [5h). We cannot rule out that this 
result is due to a still somewhat slower crossover to the 
asymptotic excluded volume dominated regime than as- 
sumed in our fit in Fig. [7j still much longer chains than 
Nf, = 50000 would be needed to check this, but this is 
a very tough task even for the PERM algorithm. On 
the other hand, we note that Eqs. (f38l) - (|43|) clearly are 
not exact, the Flory argument invariably implies that 
v = 3/5 = 0.60 instead of v = 0.588 [52| and it is un- 
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FIG. 8: Log-log plot of various possible estimates for a persistence length £ p plotted vs. qt, in both d — 2 and d = 3, as indicated. 
Here n p and (n atr ) are extracted from the use of Eq. (|71[) . cf. Fig. [5] while is taken from the fit to the initial decay of 
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amplitude £ Pt n/£b characterizing the prefactor of the asymptotic excluded volume region, {Rl) = 2£ b £ p , R N 2v with v = 3/4, 
in d — 2, versus £ p (q b )/£b- Straight line shows a fit to £ p ,nj£b = 0.3(£ p /£{,) 1//2 . (c) Log-log plot of £ p ,n/£b and N£ in d = 3 
dimensions versus £ p /£b, as extracted from the fit shown in Fig. [7] Both effective exponents and theoretical power laws are 
indicated (cf. text). 



clear to us to what extent these exponents describing the 
variation of £ Py n and N£ with t p are modified. This prob- 
lem could possibly be addressed with the renormalization 
group approach. 



V. STRETCHING SEMIFLEXIBLE POLYMERS 
IN d — 2 DIMENSIONS 

Fig. [9] presents now a selection of our results for exten- 
sion versus force curves for the two-dimensional SAW's 
of variable stiffness on a square lattice. As expected, nei- 
ther the simple result for freely jointed chains {Eq. (01} 
nor the Kratky-Porod result {Eq. (|55|) } are compatible 
with the data. For very short chains (Nb = 100) and 
intermediate values of the stiffness (qb = 0.05 which cor- 
responds to £ p ss 13£b, cf. Table HI) we note that {X)/L 
roughly agrees with Eq. (|55p . however: this agreement 
probably is not accidental, since also in the absence of 
a force for such short chains and this choice of qb the 
Kratky-Porod prediction for the mean square end-to-end 



distance (Fig. [6fb)) still is rather close to the actual re- 
sult for (Rl)/2£ p L. Similar agreement was also noted 
in our earlier work (45j | for somewhat longer and stiffer 
chains (Nb = 200 and qb = 0.03 and 0.02, respectively) 
for exactly the same reason: as long as (Rf) in the ab- 
sence of forces is still more or less correctly predicted, and 
this can be judged from the data presented in the previ- 
ous section, the general linear response relation, Eq. (0), 
which holds not only for flexible SAW's but also for stiff 
chains, implies that the K-P model still provides an ac- 
curate description of the initial linear part of the exten- 
sion versus force curve. Since in such a case where L is 
larger then £ p by only a small factor, beyond the linear 
response re gime there is no regime of Pincus blobs possi- 
ble, since y {R 2 } in Fig.QJa) and £ p then are of the same 
order (each Pincus blob needs to be formed from many 
subunits of size £ p , in order that the power law regime 
(X)/L cx (f£ p /k B T)^ 3 can develop!) Thus, we arrive 
at the general conclusion that in d = 2 the K-P result 
Eq. (|S"5"|) . is applicable only for such short chains that L 
is larger than £ p by only a small factor (L < W£ p , say), 
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FIG. 9: Relative extension (X)/L plotted versus scaled force j£ v jksT for rather flexible chains (qb = 0.4, case (a)) and for 
rather stiff chains (qt = 0.03, case (b)) in d = 2, including several different contour lengths L = Nb£t, as indicated. In (a), the 
prediction Eq. (j4}, (X) / L oc C((f£ p /kBT)(£b/£ p )), for the freely jointed chain and Eq. (|63[) for the Pincus blob prediction are 
included for comparison, while in (b) the result Eq. (|55[) for the Kratky-Porod model is included. Note that in (a) the result 
t p — 2.67ft, (Table H| was used to convert the scale from fip/ksT to f£b/ksT. Case (c) plots (X)/L versus f£ p /kBT and 
variable qt (and hence variable £ p , cf. Table U for L = 100 and case (d) for L = 400, respectively. Eq. ([55| is again included 
for comparison. 



so that in Fig. [TJa) the Pincus blob regime is essentially 
absent, and the K-P model also achieves an approximate 
description of the linear response regime. 

Of course, it is of great interest to clarify what hap- 
pens when L^> £ p . Fig. [TU] hence presents a log-log plot 
of the data for the extension versus force curves includ- 
ing long chains and rescaling the data such that a scaling 
description for the crossover from the linear response to 
the regime of Pincus blobs is obtained. One sees that 
both for flexible chains (Fig. [TOta)) and for rather stiff 
chains (Fig. |TU](b) ) a reasonable data collapse on a mas- 
ter curve is obtained, consistent with the predicted ex- 
ponents. As expected, the crossover between both power 
laws is gradual and not sharp. If one includes data for 
too large forces, one can see that the data fall system- 
atically below the Pincus power law. Similarly, when 
one puts the focus on the crossover from the Pincus blob 
regime to the saturation behavior, Fig. [TOtc). one finds 
that the data fall systematically below the Pincus power 
law for small forces (due to the crossover towards the 
linear response regime). As expected from the theoreti- 



cal considerations of Sec. II, there cannot exist a scaling 
representation which brings both crossovers of Fig. [TJa) 
to a data collapse on a master curve together. Note also 
that for the long chains the K-P model does not fit our 
data at large relative extensions (X) / L either, since our 
simulations are based on a discrete chain model. Since 
our choices of qb do not yield extremely large persistence 
lengths, the crossover from the saturation behavior pre- 
dicted by the K-P model {Eq. ([55])} to that of the FJC 
model {Eq. ([7])} is not clearly resolved either. Actually, 
for very large forces one must consider that our model is 
a lattice model, not a model of rigid bonds in the con- 
tinuum where arbitrary bond angles occurs such as the 
FJC model: hence we expect that for / — > oo the sat- 
uration behavior is 1 — (X)/L oc exp(— Jib/ksT) rather 
than k B T/f£ b . 

As a last point of this section, we consider both lon- 
gitudinal (Fig. [TTT aKb)) and transverse (Fig. [rflcHd)) 
fluctuations of the chain dimensions. These fluctua- 
tions have been normalized such that they are of or- 
der unity (and independent both L and £ p ) in the lin- 



16 




FIG. 10: Plot of y = {{X)/L)C y versus x = (jl h lk B T)C x for flexible SAW's in d = 2 (q b = 1) where the scaling factors 
C x , C v for abscissa and ordinate have been chosen C x = N 3 , C y — JVjJ , so that the coordinates (x CT ,y cr ) of the crossover 
point from the linear response regime to the Pincus blob regime in the (x,y) plane are of order unity. Several choices of L are 
included, as indicated. Part (b) is similar as (a), but for semiflexible chains with several choices of qt and L, as indicated. Now 
x = (fip/kBT)C x and the scaling factors are chosen as C x = (L/£ p ) 3 ^ 4 , C y = (L/lp) 1 ^ 4 . Straight lines in both parts indicate 
the theoretical power law (X) cx / (linear response regime) and (X) cc f 1 ' 3 (Pincus blobs regime), respectively. Part (c) is the 
same as (b), but choosing C x = Cy = 1 including only data for L = 25600, 12800, and 6400, to show the crossover from the 
Pincus blob regime to a Kratky-Porod (K-P) like regime. 



ear response regime, while in the Pincus blob regime 
(0 <C flp/ksT < 1) a crossover to a simple power law 
proportional to (/ip/fceT) 1 ' 1 ' -2 occurs. Using the same 
scaling factors C x — (L/£ p ) 3 ^ 4 (Fig. fTOT b)) for abscissa, 
a nice data collapse on the master curve is seen in Fig. [12] 
for both longitudinal and transverse fluctuations. Note 
that our scaling description for the crossover from the 
linear response regime to the Pincus blob regime, exem- 
plified in Fig. [10] and [HI does not invoke any adjustable 
parameters whatsoever (unlike the case of experiments, 
where often both L and £ p are fit parameters). However, 
the behavior at larger forces (beyond Pincus blob regime) 
is more subtle. While in the Kratky-Porod regime a 
power law proportional to (/Ip/fcsT) -2 is expected for 
large enough /, where a behavior similar as that has 
been found for freely jointed chains {Eq. (|21[) } can be 
expected, when one considers the longitudinal fluctua- 
tion {((X 2 ) - (A") 2 )/(A" 2 ) }, for the transverse fluctua- 
tion (R\) I (-Rj_)o all theories predict a slower decay (pro- 
portional to [(fip/keT)- 1 ] for large / {Eqs. (22]), pi]) . 



and (|6?p. respectively}, and this slower decay in fact is 
not seen. The reason for this discrepancy, however, prob- 
ably is the fact that in our model only kinks by ±90° are 
possible, and no small deflections are possible as in the 
K-P and FJC models. We defer a more detailed anal- 
ysis of these fluctuations to a forthcoming publication. 
Here we rather focus on the behavior of the local angu- 
lar fluctuation ((j> 2 ) (Fig. [13]) . For small forces all an- 
gles (j) between a bond and the +a:-direction are equally 
probable (at the lattice we have in d — 2 two possibili- 
ties for tf> — ir/2 or <fr = —tt/2, and two possibilities for 
(j) = or 7r, respectively). Hence for / — > we must find 
(</> 2 ) = (tt 2 +tt 2 /2)/4 = 3tt 2 /8 rj 3.7, and this is compat- 
ible with the observation. For fip/ksT > 1 we observe 
a smooth crossover towards 

{cp 2 ) cxk B T/(fe p ), fe p ^>k B T (73) 

and the crossover from ((f) 2 ) s» 3.7 to this decay seems to 
be practically independent of £ p (as shown by the super- 
position of data for different choices of qb and hence £ p 
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FIG. 11: Log-log plot of [{X 2 ) - (X) 2 ]/ (X 2 ) vs. fl p /k B T for L = 100 (a), and L = 25600 (b), including several choices for 
qt as indicated. Log-log plot of (Rj_) / {R±)o vs. /Ip/feT for L — 100 (c) and L = 25600 (d), including several choices for qt 
as indicated. A straight line with slope 1/V — 2 (y = 3/4) is shown for comparison. Data are for semiflexible chains in d = 2. 



in Fig. [LSI) and L (compare Fig. ITBT a) for L = 100 with 
Fig. HHb) for L = 25600). Although the local quan- 
tity (0 2 ) thus has a remarkably simple behavior, unlike 
the global quantities (X)/L, ((A 2 ) - (A) 2 )/(A 2 ) and 
(Rj_) I (Rj_)o, we are not aware of any theoretical predic- 
tion relating to it. 



VI. STRETCHING SEMIFLEXIBLE POLYMERS 
IN d = 3 DIMENSIONS 

We start by showing extension versus force curves 
for various choices of the contour length L in Fig. IT41 
to provide a three-dimensional counterpart to the data 
in Fig. IHl for two dimensions. It is immediately obvi- 
ous that the simple Kratky-Porod prediction {Eq. (|54p } 
does a much better job than its two-dimensional coun- 
terpart {Eq. (|55p}. Again, we emphasize that there are 
no adjustable parameters whatsoever in our comparison, 
L = Nfjlb is trivially known, and £ v comes from Fig. EJd) . 
In fact, for L = 200 and L = 400 most of the data for 
0.01 < (A) < 0.3 follow the K-P prediction, for a wide 
range of choices for g& and hence l v {Table ITT]}- . only data 
for rather flexible chains (such as qb — 0.4, for which 
£ p = 1.2) deviate strongly from the K-P model, as ex- 



pected. The simple Langevin function {Eq. (HJ)} does 
not describe the behavior of these lattice chains with dis- 
crete bond angles. Note that these chains are too short 
to show a well-developed Pincus blob regime yet. For 
(A) /L > 0.3 systematic deviations from the K-P predic- 
tion occur, which we attribute to effects due to the dis- 
creteness of bonds and bond angles in our model. Only 
for very long chains (such as N = 6400 and 25600) do we 
find more pronounced deviations from the K-P prediction 
also for small relative extensions, (A) / L < 0.1. 

For long chains, however, we do expect to see excluded 
volume effects (manifested in Pincus blobs) , as discussed 
in Sec. II. Thus Fig. [13] presents our data in suitably 
scaled form, considering the crossover from the linear re- 
sponse regime to the Pincus blob regime, both for flex- 
ible chains (a) and semiflexible ones (b), as well as the 
crossover from Pincus blobs to Kratky Porod behavior 

(c) , and we show a close-up of the Kratky-Porod regime 
for very long and at the same time rather stiff chains 

(d) . Of course, even with chain lengths up to 25600 
it is not yet possible to clearly resolve all the different 
power laws shown in Fig. [1] (b): in order to be able 
to distinguish the various crossovers clearly from each 
other, we would need very stiff chains (£ p should then be 
in the range 10 2 < £ p /£b < 10 4 ), and then one would 
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FIG. 12: Log-log plot of [(X 2 ) - (X) 2 ]/(X 2 ) (a) and {R\)/{R\)o (b) vs. (fl p /k B T)C x with C x = (L/£ p ) 3/4 (see Fig^b)) 
for semiflexible chains in d — 2 including several choices of L and qb, as indicated. 



(a) 



(b) 





fL/I^T 



FIG. 13: Log-log plot of (<j> 2 } vs. fip/ksT for L — 100 (a), and L = 25600 (b), including several choices for qt, as indicated. 
((j> 2 ) = 3tt 2 /8 as / ->■ 0. Data are for semiflexible chains in d = 2. 



need to have chain lengths of many millions in order 
to have a well-developed Pincus-blob regime. Thus, we 
can verify the Pincus blob regime only for rather flexi- 
ble chains (Fig. [TBT a)). for which then a well-developed 
Kratky-Porod regime is absent. For the stiff chains, we 
can see some tendency of the data to deviate from the 
K-P regime in the direction towards the Pincus blob 
regime (Fig. ll5f b)(c)). but the latter is not fully reached 
because the crossover to the linear response takes over 
(Fig. ll5f b)). And when we study very stiff chains, we find 
deviations from the K-P model for rather small (X)/L 
already, due to the discrete character of our chains. 

Fig. [T5] shows again data for the normalized fluctua- 
tions of the chain linear-dimensions, and Fig. 1 171 presents 
a counterpart to Fig. [T5J showing a log-log plot of the 
local fluctuation (<fi 2 ) versus fip/ksT. While the latter 
(for large L) show again a simple crossover from the con- 
stant 7r 2 /3 describing (</> 2 ) for small forces to a power law 
ksT/ftp for f£p/ksT > 1, as in d = 2 dimensions, the 
behavior of the fluctuations in the chain linear dimen- 
sions clearly is rather complicated. Of course, there is a 



need to extend the scaling analysis, that was presented 
for (X) as a function of fip/ksT in Fig.QJb) to the fluc- 
tuations (X 2 ) — (X) 2 and (Rj_) in greater detail than we 
have done so far. We expect that analyzing these fluctua- 
tions should yield additional and valuable information on 
the structure of stretched semiflexible chains, and allow 
to pin down the parameters needed to relate experimen- 
tal data to theoretical models more precisely. We plan 
to tackle this task in a forthcoming study. 



VII. CONCLUSION 

In this paper, we have studied self- avoiding walks on 
square and simple cubic lattices, where an energy penalty 

associated with chain bending to model semiflexibility 
of the polymer chains, by extensive Monte Carlo simu- 
lations, using the PERM algorithm. We have obtained 
both force versus extension curves and chain linear di- 
mensions in the absence of forces for a wide range of chain 
lengths Nh (typically iVf, up to 25600) and chain stiff- 
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FIG. 14: Log-log plot of (X)/L versus flp/ksT for several choices of q b as indicated, for contour length L = 200 (a), L = 400 
(b), L = 6400 (c), and L = 25600 (d). Full curve always refers to the K-P model prediction, Eq. (|54[) . Broken curve in (a) is 
the Langevin function, Eq. (f4]), using t v /l b = 0.71 for q b = 0.4 (Table Ell- Data are for semiflexible chains in d = 3. 



ness (characterized by qj, = exp(— Cb/ksT)). In Sec. II, 
we have attempted to present a coherent phenomcnolog- 
ical theoretical description, combining results from scal- 
ing concepts with other results derived from the Kratky- 
Porod model, to explain the various crossovers that can 
occur in the force versus extension curves for various cir- 
cumstances (Fig.Q}. We have emphasized that the case of 
d = 2 dimensions is rather different from the case d = 3: 
only in the latter case one can identify a linear regime in 
the force versus extension curve that is compatible with 
the Kratky-Porod model; the linear response regime both 
in d = 2 and d = 3 dimensions is strongly affected by the 
presence of excluded volume effects, and for long enough 
chains is followed by a nonlinear ( "Pincus blob" ) regime 
for stronger forces both in d = 2 and d = 3 dimensions. 
However, for very stiff and not too long chains in d = 3 
the chains in the absence of a force show Gaussian behav- 
ior, and in this case the stretched chains do not exhibit 
the nonlinear Pincus blob regime, and the Kratky-Porod 
model holds throughout (apart from very strong forces, 
where the discrete character of polymer chains matter). 

The Monte Carlo data that we have generated do pro- 
vide evidence for these concepts, particularly in the rel- 
atively simple case of d = 2 dimensions. While in d = 2 
all expected regimes of the force vs. extension curves are 



confirmed and the expected scaling behavior is verified, 
problems remain concerning the precise understanding 
of longitudinal and transverse fluctuations of chain lin- 
ear dimensions of the chains. More work on these as- 
pects (from theory, simulation, and experiment) clearly 
is desirable. We recall that imaging techniques can pro- 
vide rather detailed information on chain configurations 
of semiflexible polymers adsorbed on substrates; we ex- 
pect that our work should be useful to interpret such 
experiments. 

In the case of d = 3 dimensions, our numerical evidence 
is much more limited: chain lengths Nb — 25600 clearly 
do not suffice to fully resolve three distinct power laws 
(separated by smooth crossovers) in the force versus ex- 
tension curves. However, simulations for chains that are 
one or two orders of magnitude larger clearly are not fea- 
sible at present. We do obtain evidence, however, that 
for rather stiff thin short chains excluded volume effects 
indeed are negligible, as expected, and hence the Kratky- 
Porod model holds. However, when chain stiffness is due 
to thickness (persistence length i p being proportional to 
local chain diameter D), the regime of Gaussian statistics 
disappears and rather excluded volume effects dominate 
throughout, resulting in a rather broad regime where the 
force versus extension curve is nonlinear already for small 
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FIG. 15: (a) Log-log plot of (X) / (£ p £ b L 3 ) 1/5 vs f/f c = fL 3/5 {£ p l b ) 1/5 /k B T, where f c is the crossover force defined in Eq. (60}, 
remembering D = lb in our model, for the choices q b = 1.0, 0.4, and 0.2 and many choices of L (a) and for q b = 0.1, 0.05, 
0.03 but only L = 25600 and L = 12800 (b). Case (c) shows a plot of ((X) /L)(£ p /£ b ) versus (f£ p /k B T)(£ p /£ b ) for the choices 
q b = 0.05, 0.1, and 0.2, again for L — 25600 and L — 12800 only, to test for indications of a crossover from Pincus blobs to 
the Kratky-Porod model, showing only the vicinity of the region when this crossover should occur. Case (d) is a blow up the 
region 0.05 < (X)/L < 0.5, for q b = 0.02, 0.01, 0.05, L = 25600 and L = 12800, to show the full K-P region for rather long and 
rather stiff chains. Data are for semiflexible chains in d — 3. 



(X)/L. Also for rather flexible chains, clear evidence for 
the Pincus blob regime is obtained (Fig. [rSla 1 ) 1 ). 

In our modelling, we have approximated the interac- 
tions between monomers of the chain as a strictly lo- 
cal excluded volume interaction. Of course, in many 
cases of interest the interactions are of longer range, 
e.g. because of electrostatic interactions between charged 
groups. Particularly for polyelectrolytes the resulting 
problem of an "electrostatic persistence length" has re- 
ceived longstanding attention in the literature 
[79j . Molecules such as DNA and RNA do possess a sub- 
stantial linear charge density, and the properties of such 
polyelectrolytes, in fact, will depend on the electrostatic 
screening due to ions in the solution, and thus the effec- 
tive persistence length will depend on ionic strength. In 
this context, the concept of an "effective thickness" of 
polyelectrolytes, that are described in terms of a "thick 
chain model" [79|], has been used to model experimental 
extension versus force curves. It will be an interesting 
task for the future, beyond the scope of the present pa- 
per, to clarify the extent to which such approaches are 



equivalent to the scaling concepts applied here. In any 
case, it is very reassuring that very recently, after our 
study was completed, single-molecule elasticity measure- 
ments of the onset of excluded volume effects of stretched 
poly(ethylene glycol) were published 80]. In this work, 
the Pincus blob scaling behavior (L cx f 2 ^ 3 in d = 3 
dimensions) could be seen under several circumstances, 
followed by a crossover to the linear behavior (L oc /) 
and subsequent saturation, compatible with the behav- 
ior predicted by the K-P model. The interpretation given 
in Ref. [80] for these experiments is fully consistent with 
the description given in the present paper. Some ear- 
lier evidence for the Pincus behavior was also found for 
DNA HI, HI]. 

In our simulations, we have not considered the ef- 
fects of varying the diameter D of our chains (we have 
chosen D = £{, = 1, the lattice spacing, throughout). 
For biopolymers D is a parameter of great interest as 
well 1821. of course. In our studies of bottlebrush poly- 
mers ^3,13], however, we studied conditions for which 
D oc £ p , and then the K-P model was not useful even in 
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FIG. 16: Log-log plot of [(X 2 ) - (X} 2 ]/(X 2 } vs. f£ p /k B T for L = 100 (a) and L = 25600 (b), and for various choices of q b , 
as indicated. Part (c) shows analogous data for (i? 2 L )/(i? 2 )o vs. fip/ksT for L = 100 and part (d) for L — 25600. A straight 
line with slope 1/u — 2 (y = 0.588) is shown for comparison. Data are for semifiexible chains in d — 3. 
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FIG. 17: Log-log plot of ((f> 2 ) vs. fi p /k B T for L = 100 (a) and L = 25600 (b), and for many choices of qt, as indicated. 
((f) 2 ) = 7r 2 /3 as / — > co. Data are for semifiexible chains in d — 3. 



d = 3 dimensions. 

Thus we hope that the present work will contribute to 
the better understanding of both existing and future ex- 
periments. A very interesting aspect, completely beyond 
the scope of the present work, are dynamic properties 
of stretched semifiexible polymers in solution, see e.g. 
Ref. [83]. Our study should yield useful inputs for such 



problems, too. 
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